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Abstract 

The scaling behaviour of the persistence probabihty in the critical dy- 
namics is investigated with both the heat-bath and the Metropolis algo- 
rithm for the two-dimensional Ising model and Potts model. Special atten- 
tion is drawn to the dependence on the initial magnetization. The global 
persistence exponent is measured. Universality is confirmed. 
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For long a set of two static exponents /3, v and one dynamic exponent z 
has been used to characterize the scahng behaviour of critical dynamics. This 
is sufficient when the dynamic systems almost reach the equilibrium. When 
the dynamic systems are far from the equilibrium, it was believed that the be- 
haviour of the dynamic systems essentially depends on the microscopic details 
and there exists no universal scaling behaviour. However, it has recently been 
discovered that for a dynamic process, in which a magnetic system initially at 
very high temperature is suddenly quenched to the critical temperature and then 
released to the dynamic evolution of model A, already in the macroscopic short- 
time regime emerges universal scaling behaviour. A new independent exponent 
should be introduced to describe the dependence of the scaling behaviour on the 
initial magnetization [0, ^, §]. Furthermore, if a non-zero but small magnetiza- 
tion is given to the initial state, at the macroscopic early time the magnetization 
surprisingly undergoes a critical initial increase [|I], H, ^ |], ^ H 

M(^)~mot^ (1) 

where 6 is the new critical exponent which is related to the dimension xq of the 
initial magnetization mo by ^ = {xq — (3/v)/ z. 

Are there any more critical exponents in the dynamic system described above? 
Very recently, it was argued that the global persistence exponent of the magne- 
tization is another new independent critical exponent p. This relies on the 
observation that the dynamic evolution of the magnetization is not a Markovian 
process. If a non-zero magnetization is given to the initial state at very high 
temperature, the time-dependent magnetization keeps its sign at the time t with 
a probability p{t). In the limit of zero initial magnetization, the probability p{t) 
decays by a power law with respect to the time t 

p{t)-t-'\ (2) 

Here 9i is called the global persistence exponent. The exponent 9i has been 
calculated perturbatively for the 0{N) vector model and numerically determined 
for the two-dimensional Ising model with the heat-bath algorithm 0. In reference 
0, the power law behaviour in Eq. (0) was not clearly observed since the lattice 
size was not big enough. Only a combination of the exponents diz has been 
obtained from finite size scaling. The result carried a relatively big error. Stauffer 
has performed simulations with the heat-bath algorithm for bigger lattice sizes 
and estimated the exponent 6i from the power law behaviour for the two- and 



three-dimensional Ising model [0. For the two-dimensional Ising model with 
lattice size L = 300, the result is 6i = 0.225. The error is probably somewhat 
below 0.01 |Ty]. If the time evolution of the magnetization is a Markovian process, 
a scaling law which relates both the exponent 6 and 6i 

e^ = ai = -e + {l- (3/u)/z (3) 
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can be deduced. It has been argued at the two-loop level of an e-expansion that 
this scaling law is violated . It is interesting to investigate numerically how big 
the violation of the scaling law is. Due to the relative big errors for the exponent 
01 obtained in references ^, the results are not completely satisfactory]^. 

In this letter, we present systematic Monte Carlo simulations for the scaling 
behaviour of the persistence probability of the magnetization. Special attention 
is drawn to the dependence of the scaling behaviour on the initial magnetization. 
Calculations are extended to the two-dimensional Potts model. Further, to con- 
firm universality, all the simulations are carried out with both the heat-bath and 
Metropolis algorithm. 

At first, let us concentrate on the two-dimensional Ising model. In reference 
@, the initial state is simply generated randomly — with equal probability of 
spins up and down. Since the lattice size in practical simulations is finite, it 
is almost always the case that a finite non-zero initial magnetization will be 
generated in each initial configuration. Then the system is released to evolve at 
the critical temperature with the heat-bath algorithm. We update the system 
until the magnetization change its sign. Repeating this procedure many times 
with different initial configurations and random numbers, we can measure the 
global persistence probability p(t). If the lattice is sufficient large, p{t) should 
decay by the power law given in Eq.(0). As was observed in reference Q, however, 
the lattice size L = 128 is not sufficiently big and the power law is not clearly 
observed due to the finite size effect. 

Here we should point out that in this dynamic system there are two kinds of 
finite size effects. One is the normal finite size effect which takes place in a time 
scale ti ~ L^- Whenever the system evolves into this time regime, p{t) will decay 
by an exponential law exp(— t/t^.) rather than a power law. For a lattice size as 
L = 128 or bigger, the time scale ti is quite big. Later we will see that in the 
short-time regime this effect is not prominent. Another kind of finite size effect 
is an extra finite size effect from the initial configurations. When the lattice 
size is finite, each initial configuration generated randomly as described above 
has a finite non-zero initial magnetization. This non-zero initial magnetization 
tends to zero only in the limit of infinite lattice size. Such a finite non-zero 
initial magnetization would modify the power law behaviour in Eq.(0) even in 
the short-time regime of the dynamic evolution. 

In order to handle the extra finite size effect efficiently and to discuss the 
dependence of the scaling behaviour on the initial magnetization, we introduce 
a sharp preparation technique in preparing the initial configurations: first the 
initial configuration is generated randomly; then we select a site randomly on the 
lattice and flip the spin if the updated magnetization comes nearer to a pre-fixed 
small value mo', we continue this procedure until the initial magnetization reaches 

^When the revised version of this paper was completed, the two loop contribution has been 
explicitly calculated 
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the value niQ. By this technique, we can obtain any fixed value for the initial 
magnetization. 

In Fig. |l], the persistence probability p{t) is displayed with solid lines in double- 
log scale for the lattice size L = 128 and different initial magnetization rriQ. For 
comparison, the result without the sharp preparation of the initial magnetization 
is also plotted with a dotted line. Samples of independent initial configurations 
are 20 000, 40 000 and 60 000 for mo = 0.01, 0.001 and 0.0005 respectively and 
that without the sharp preparation of the initial magnetization is 40 000. From 
the figure, we see clearly that the behaviour of p{t) apparently depends on the 
initial magnetization. When the initial magnetization tends to zero, however, 
it convergences. The difference between the curves with mo = 0.001 and mo = 
0.0005 is already very small. The curve without the sharp preparation technique 
is a certain combination of those with different mo. For the lattice size L = 128, 
the curves with small mo do also not yet present perfect power law behaviour 
even though it is improved compared with that without the sharp preparation 
technique. 

In Fig. 0, the persistence probability p{t) is displayed in double-log scale for 
different lattice sizes. The sharp preparation technique is adopted, mo is taken 
to be 0.0005. The dashed line, upper solid fine and the dotted line correspond to 
the lattice size L = 128, 256 and 512 respectively. Samples of the independent 
initial configurations are 60 000 and 28 000 for the lattice sizes L = 256 and 512. 
In the figure, it can be seen that after a microscopic time scale tmic ~ 50, the 
power law decay is observed for bigger lattices. If we measure the slopes of these 
three curves from a time interval [50, 1000], the resulting exponent 6i is 0.224(2), 
0.238(3) and 0.238(6) respectively. The errors are estimated by dividing the data 
into three or four groups. As pointed out above, there is still some extra finite 
size effect for the lattice size L = 128, however, the lattice size L = 256 and 
L = 512 give the same results within the statistical errors. The extra finite size 
effect is already negligible small. 

In Fig. 0, we can also see that the normal finite size effect has not yet appeared 
up to the evolution time t = 1000. The normal finite size effect is characterized by 
an exponential decay of p{t), which is faster than the power law decay. However, 
for the smaller lattice size, e.g. L = 128, the practically measured exponent is 
even smaller. Such a tendency can also be seen clearly in the figure. 

Are the power law behaviour of p{t) and the critical exponent 6i universal? In 
order to check this point, we have also performed simulations with the Metropolis 
algorithm. For smaller lattice sizes, the behaviour of p{t) with the heat-bath and 
the Metropolis algorithm is somehow different. When the lattice is sufficient 
large, however, both algorithms present almost the same power law behaviour. 
The result of p{t) for the Metropolis algorithm with lattice size L = 256 and 
mo = 0.0005 is displayed with the lower solid line in Fig. ^ The corresponding 
exponent is 9i = 0.236(3). Within the statistical errors, the result is consistent 
with that of the heat-bath algorithm. Universality is confirmed. 
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To have more confidence on our measurements, we have also performed simu- 
lations with the heat-bath algorithm for the lattice size L = 256 with mo = 0.001 
and with initial magnetization not sharply prepared. The results are displayed 
with a dashed line and a solid line in Fig. |^. The measured exponents are 0.237(5) 
and 0.233(5) respectively. The result of mo = 0.001 is the same as that of 
mo = 0.0005 within the errors. The result without the sharp preparation tech- 
nique is slightly smaller, but the extra finite size effect is for L=256 already not 
so big. Our result 9i = 0.233(5) without the sharply prepared initial magneti- 
zation, i.e. with a random initial state, is slightly different from what Stauffer 
obtained, di = 0.225(10). But the difference is still within the stattistical errors. 
For comparison, the result for L = 256 and mo = 0.0005 is also plotted with the 
dotted line in Fig. ^. 

All the results for the Ising model have been summarized in Table ^ and 
Table |[ Samples of the independent initial configurations are from 28 000 to 
60 000, depending on the lattice size L and the initial magnetization mo. Errors 
are estimated by dividing the samples into three or four groups. Compared with 
the result obtained from the finite size scaling in reference 0, our results are 
more accurate. 

Encouraged by the success for the Ising model, we have carried out similar 
simulations for the two-dimensional Potts model. In Fig. ^, p{t) has been plotted 
in double-log scale with mo = 0.0005. The dashed line and the upper solid line 
are the result with the heat-bath algorithm for lattices L = 144 and 288, while 
the dotted line and the lower solid line are those with the Metropolis algorithm. 
Slightly different from the case of the Ising model, for the lattice size L = 144 
the extra finite size effect is already not so big. For example, with the heat-bath 
algorithm the measured exponent 6i is 0.344(6) and 0.350(1) for the lattice sizes 
L = 144 and L = 288 respectively. Within the statistical errors, they cover each 
other. In Table 0, the exponent 6i measured for the Potts model for the lattice 
size L = 288 with both the heat-bath and the Metropolis algorithms is also given. 
The statistics is 60 000. 

How big is the violation of the scaling law in Eq.(^) in the numerical sense? 
The static critical exponent (3 and u for the two-dimensional Ising and Potts 
model are exactly known, while the dynamic exponent 6 and z have numerically 
been measured from the short-time dynamics 0, ^. The value of ai = —9 + (1 — 
l3/v)/z for the Ising model and the Potts model with both the heat-bath and the 
Metropolis algorithm are listed in Table ^ These results are taken from reference 
[[7|. Apparently the exponent 6i is not equal to the exponent ai. The difference is 
about 10 percent. For example, for the Ising model with the heat-bath algorithm 
Oi = 0.238(3), while ai = 0.215(1). 

Here we should mention that in the estimate of the exponent ai both ex- 
ponents 6 and z are involved. The numerical measurements of 6 are somehow 
satisfactory 0, The value of 6 for the Potts model is also relatively small and 
will not induce a big error for ai. However, the exponent z is problematic, at 
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least for the Ising model. Values ranging from z = 2.13 to z = 2.17 can be found 
in the literature |jT2|, |1^, |1^, |15|, |16|, |T^. Fortunately, the value of ai is not so 



sensitive to the exponent z. For example, in reference we find 9 = 0.191(3) 
and z = 2.172(6) for the Ising model with the heat-bath algorithm. Then one 
obtains ai = 0.212(3). It is slightly smaller than that given in Table H, but still 
consistent within the statistical errors. 

In reference 0, a combination of the exponents 6iz = 0.505(20) was obtained 
from finite size scaling. This leads to the values 6i = 0.233(9) and 6i = 0.234(9) 
for z = 2.172(6) and z = 2.155(3) respectively. These values are smaller compared 
with those in Table |^, and taking into account the bigger error, the difference 
between the exponent 6i obtained here and ai given in Table Q is less prominent. 
- In reference |]lOl the result 6i = 0.225 with an error around 0.01 is given. Here 



the values for 6i and ai nearly cover each other. 

Finally, we investigate the general scaling behaviour for the persistence prob- 
ability p{t) in case the effect of the initial magnetization can not be neglected. 
Let us assume a finite size scaling form 

p{t,L,mo) =r'^f{t'/^L-\e<^/^mo). (4) 

with xo being the scaling dimension of the initial magnetization mo. From the 
time evolution of the magnetization or the auto-correlation, the value of xq for the 
Ising and Potts model has previously been determined directly or indirectly 
through its relation with the exponent 9, i.e. Xq = 6^2; + /5/z/ P, |^. 

Now, we intend to determine the scaling dimension xq from the scaling form 
in Eq. (^. In order to do this, we take the two-dimensional Ising model as an 
example and perform a simulation for the lattice size Li = 256 with the initial 
magnetization moi = 0.0050. If the scaling form in Eq. (^ holds, one should be 
able to find an initial magnetization mo2 for the lattice size L2 = 128 such that 
the curves of p{t) for both lattice sizes collaps. From the ratio of moi and mo2, 
one may estimate Xq. Practically we have performed simulations for L2 = 128 
with two values of the initial magnetization, mo = 0.0090 and mo = 0.0115. By 
linear extrapolation, we can obtain data with the initial magnetization between 
these two values of mo. We search a curve which has a best fit to the curve for 
Li = 256 and then determine mo2. In Fig. |^, such a scaling plot is displayed. 
The lower and upper solid lines are the persistence probability for L2 = 128 
with mo = 0.0090 and mo = 0.0115 respectively, while the crosses represent the 
rescaled one for Li = 256 with moi = 0.0050. The solid line fitted to the crosses 
is the persistence probability for L2 = 128 with mo2 = 0.0101(1). From the 
figure, we see clearly the data collaps nicely and we can conclude that the scaling 
form in Eq. (^) holds. The estimated scaling dimension for mo is Xq = 1.01(1). 
To confirm this result, we have also done similar simulations for Li = 256 with 
moi = 0.01 and obtained Xq = 1.06(4). Both values are consistent and the finite 
mo effect is already quite small. Further, the simulations with the Metropolis 
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algorithm show similar results. However, the value of xq measured here is very 
different from that measured from the time evolution of the magnetization or the 
auto-correlation, xq = 0.536(2) 0. This fact remains to be understood. 

In conclusion, we have numerically measured the global persistence probability 
p{t). The critical exponent 6i is directly determined for the two-dimensional Ising 
and Potts model with both the heat-bath and the Metropolis algorithm. Within 
the statistical errors, universality is confirmed. Compared with the previous 
results for the Ising model with the heat-bath algorithm, our simulations have 
been carried out with special attention to the dependence of the scaling behaviour 
on the initial magnetization. 
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L 


256 


mo 


random 


.0010 


.0005 


9i 


.233(5) 


.237(5) 


.238(3) 



L 


512 


mo 


.0005 


01 


.238(6) 



Table 1: The exponent 6i measured for the two-dimensional Ising model with 
the heat-bath algorithm. 'Random' represents the result obtained without the 
sharp preparation technique. 





Ising 


Potts 




HB 


ME 


HB 


ME 


ei 


.238(3) 


.236(3) 


.350(1) 


.350(8) 




.215(1) 


.212(2) 


.320(3) 


.324(3) 



Table 2: The exponent di measured for the two-dimensional Ising and Potts 
model with both the heat-bath (HB) and the Metropolis (ME) algorithm. The 
initial magnetization is mo = 0.0005 and the lattice sizes are L = 256 and 288 
for the Ising and Potts model respectively. The exponent ai — —9 + {1 — l3/v)/ z 
is taken from reference [7]. 
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Figure 1: The persistence probability p{t) for the Ising model with the heat- 
bath algorithm. The lattice size is L = 128. Sohd lines are the results with the 
sharp preparation of the initial magnetization. The initial magnetization mo is 
0.01, 0.001 and 0.0005 (from above). The dotted line is that without the sharp 
preparation technique. 
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Figure 2: The persistence probability p{t) for the Ising model with the heat- 
bath algorithm and the Metropolis algorithm. The sharp preparation technique 
is adopted with mo = 0.0005. The dashed line, upper solid line and the dotted 
line are the results of the heat-bath algorithm for the lattice sizes L = 128, 256 
and 512 respectively. The lower solid line is that of the Metropolis algorithm for 
the lattice size L = 256. 
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Figure 3: The persistence probability p{t) for the Ising model with the heat-bath 
algorithm. The lattice size is L = 256. The dotted line and the dashed line are 
the results for mo = 0.0005 and 0.001 respectively, while the solid line shows the 
result for the initial magnetization not sharply prepared. 
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Figure 4: The persistence probability p{t) for the Potts model with the heat- 
bath algorithm and the Metropolis algorithm. The sharp preparation technique is 
adopted and itlq = 0.0005. The dashed line and the upper solid line are the results 
of the heat-bath algorithm for the lattice sizes L = 144 and L = 288 respectively. 
The dotted line and the lower sohd line are the results of the Metropolis algorithm 
for the lattice sizes L — 144 and L — 288. 
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Figure 5: The scaling plot of the persistence probability p{t) for the two- 
dimensional Ising model with the heat-bath algorithm. The lower and upper 
solid lines are the data for L2 = 128 with mo = 0.0090 and mo = 0.0115 respec- 
tively. The crosses are those for Li = 256 with moi = 0.0050. The solid line 
fitted to the crosses corresponds to L2 = 128 and mo2 = 0.0101(1). 
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